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Abstract — High  frequency  (HF)  interference  in  surface  wave 
over-the-horizon  (OTH)  radars  typically  exhibits  a  time-varying 
or  non-stationary  spatial  structure.  Adaptive  beamformers  that 
update  the  spatial  filtering  weight  vector  within  the  coherent 
processing  interval  (CPI)  are  likely  to  suppress  such  interference 
most  effectively,  but  the  intra-CPI  antenna  pattern  fluctuations 
result  in  temporal  de-correlation  of  the  clutter  which  severely 
degrades  sub-clutter  visibility  after  Doppler  processing.  A  robust 
adaptive  beamformer  that  effectively  suppresses  non-stationary 
interference  without  degrading  sub-clutter  visibility  is  proposed. 
The  new  algorithm  is  computationally  efficient  and  suitable 
for  practical  implementation.  Its  operational  performance  is 
evaluated  using  experimental  data  recorded  by  the  Uuka  HF 
surface  wave  (HFSW)  OTH  radar,  located  near  Darwin  in  far 
north  Australia. 

1.  Introduction 

HFSW  radars  operate  by  normal  line-of-sight  propagation 
in  addition  to  surface  wave  propagation  [1],  The  latter  enables 
surface  and  airborne  targets  to  be  detected  at  ranges  of  up  to 
200-300  kilometres  over  sea  water  [2],  [3].  Offshore  targets 
of  emerging  interest  include  “sea-skimming”  missiles  [4],  and 
small  “go-fast”  boats  used  for  smuggling  and  drug  trafficking 
[5],  [6].  Doppler  processing  can  usually  resolve  high  velocity 
targets  from  the  most  powerful  (first-order)  ocean  clutter 
returns.  This  often  allows  fast-moving  targets  with  small 
radar  cross  section  (RCS)  to  compete  for  detection  against 
interference-plus-noise. 

In  the  user-congested  HF  spectrum,  unoccupied  frequency 
channels  of  sufficient  bandwidth  (50-100  kHz)  are  at  times 
extremely  difficult  to  find,  particularly  in  the  lower  HF  band 
(5-15  MHz)  or  at  night  time  when  the  ionosphere  is  known 
to  propagate  interference  sources  very  long  distances.  The 
ionospheric  or  skywave  propagation  is  detrimental  to  HFSW 
radar  because  it  effectively  increases  the  number  of  interferers 
at  the  radar  site,  making  it  sometimes  impossible  to  find  a 
“silent”  frequency  channel  to  operate  the  radar  [7],  [8].  In  such 
situations,  the  signal-to-interference  plus  noise  ratio  (SINR) 
can  be  enhanced  by  using  adaptive  beamforming  to  cancel 
directional  interference  from  natural  and  man-made  sources. 

Radio  frequency  interference  (RF1)  typically  exhibits  a  time- 
varying  spatial  structure  statistically  characterised  by  a  non- 
stationary  spatial  covariance  matrix  over  the  CPI.  This  physical 


phenomenon  may  arise  due  to  the  dynamic  properties  of  the 
ionospheric  layer(s)  propagating  interference  [9],  the  variation 
in  geometry  between  radar  receiver  and  interference  source(s) 
[10],  or  the  impulsive  nature  of  the  source(s)  [11].  Experi¬ 
mental  investigations  confirm  that  RF1  spatial  non-stationarity 
within  typical  OTH  radar  CPI  has  a  pronounced  impact  on 
adaptive  beamformer  performance  and  must  be  taken  into 
accounted  to  ensure  effective  HF  interference  suppression,  see 
[9],  1 12]  for  example.  The  long  dwell  times  required  by  HFSW 
radars  (CPIs  of  30-120  seconds)  means  that  standard  adaptive 
beamforming  implementations  based  on  time-invariant  array 
weight  vectors  often  fail  to  significantly  improve  the  SINR. 
For  this  reason,  such  systems  are  highly  susceptible  to  the 
RFI  spatial  non-stationarity  phenomenon. 

The  alternative  is  to  update  the  adaptive  beamforming 
weights  within  the  CPI.  Time-varying  weights  are  likely  to 
yield  more  effective  interference  cancellation  while  providing 
distortionless  (unity  gain)  response  to  target  signals  incident 
from  the  radar  look  direction,  but  the  beamformed  clutter 
component  may  suffer  a  severe  loss  of  temporal  coherence  as 
it  is  received  over  the  entire  adaptive  antenna  pattern  which 
(for  non-stationary  interference)  may  be  required  to  fluctuate 
during  the  CPI.  This  degradation  manifests  itself  as  a  smearing 
of  the  clutter  energy  across  the  Doppler  (velocity)  search  space 
with  the  potential  to  obscure  the  presence  of  useful  signals. 

We  propose  a  robust  adaptive  beamformer  that  effectively 
suppresses  non-stationary  interference  without  significantly 
degrading  sub-clutter  visibility.  A  desirable  property  of  the 
algorithm  is  that  its  computational  cost  is  much  lower  than  the 
Stochastic  Constraints  (SC)  method  previously  used  to  address 
this  problem  [14].  Apart  from  its  computational  efficiency,  the 
new  method  is  also  more  robust  to  estimation  errors  that  can 
accumulate  with  each  update  and  degrade  Doppler  processing 
performance.  Both  qualities  make  it  an  attractive  candidate  for 
practical  implementation  in  operational  HFSW  radar  systems. 
The  array  data  model  is  described  in  section  2.  Section  3 
describes  the  proposed  algorithm.  Section  4  presents  and  dis¬ 
cusses  the  experimental  results  that  compare  the  performance 
of  the  proposed  algorithm  with  the  conventional  beamformer, 
the  sample  matrix  inverse  (SMI)  technique  and  the  SC  method. 
Conclusions  are  given  in  section  5. 
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II.  Data  Model 

Let  us  define  x(f)  as  the  complex  iV-dimensional  antenna 
array  snapshot  vector  received  at  a  particular  range  cell  in 
a  CPI  composed  off  =  0,l,...,P  —  1  pulses.  In  general, 
x(f)  contains  an  additive  mixture  of  radar  clutter  c(f),  external 
interference  i(t),  internal  noise  n(f)  and  potentially  a  target 
signal  s (f): 

x(f)  =  s(f)  +  c(f)  +  i(f)  +  n(f)  (1) 

A  far-held  target  echo  received  by  a  narrowband  uniform  linear 
array  (ULA)  can  be  expressed  as  s(t)  =  g(t)v(6 ),  where  the 
scalar  waveform  g(t)  represents  the  target’s  temporal  signature 
and  v(0)  is  the  spatial  response  (steering  vector)  for  an 
angle-of-arrival  9  that  coincides  with  the  radar  look  direction. 
An  ideal  target  of  constant  reflectivity  and  radial  velocity 
throughout  the  CPI  has  a  signature  g(t)  =  where 

jie^  is  the  complex  amplitude  and  fr[  is  the  Doppler  frequency 
shift  normalised  by  the  pulse  repetition  frequency  (PRF). 
Under  these  conditions,  the  simplest  target  signal  model  may 
be  defined  as 

s(t)  =  vej^fdt+*\l,ej2n$sin0,...,e2*(N-1)*sin0]T  (2) 

where  d  is  the  spacing  between  adjacent  antenna  elements, 
A  is  the  operating  wavelength  and  (-)T  denotes  transpose. 
The  internal  receiver  noise  is  assumed  to  be  uncorrelated  with 
itself  over  different  receivers  and  pulses.  It  is  spatially  and 
temporally  white  with  the  following  correlation  properties: 

E{n(r)nff(s)}  =  Srs  (3) 

where  of  is  the  noise  power,  Srs  is  the  Kronecker  delta,  E{-} 
is  the  statistical  expectation  and  H  denotes  the  Hermitian 
transpose. 

A.  External  Interference 

Let  R(f)  =  E{(i(f)  +  n(t))(i(t)  +  n(f))H}  be  the 
interference-plus-noise  spatial  covariance  matrix  at  time  t, 
the  variation  across  different  pulse  repetition  intervals  (PRIs) 
arises  only  from  the  spatial  non-stationarity  of  the  external 
interference  i(t),  that  is, 

R(f)  =  E{i(f)iH(f)}  +  a2I  (4) 

As  R(f)  is  unknown  in  practice,  it  is  often  replaced  by 
the  sample  spatial  covariance  matrix  R(f)  estimated  using 
secondary  (clutter-free)  snapshots  jfc(f)  =  i kit)  +  n &(£). 
Usually,  secondary  data  is  extracted  from  a  limited  number 
of  K  clutter-free  range  cells  available  at  the  tth  pulse.  Using 
these  snapshots,  the  sample  estimate  of  the  interference-plus- 
noise  covariance  matrix  is  given  by 

R-(*)  =  {t)  (5) 

fc= l 

In  a  bistatic  HFSW  radar,  secondary  data  may  be  extracted 
from  range  cells  prior  to  the  reception  of  the  direct  wave.  In 
the  monostatic  case,  the  high  attenuation  of  the  surface  wave 
with  distance  (more  severe  than  the  inverse  square  law)  allows 


the  clutter-free  cells  to  be  obtained  at  sufficiently  long  range. 
Otherwise  the  clutter  must  be  pre-filtered  [15],  [16],  [17]. 

If  the  CPI  is  partitioned  into  M  sub-CPI  or  data  batches 
containing  Q  pulses  (M  =  P/Q),  the  adaptive  beamforming 
weight  vectors  wm  for  to  =  1  can  be  updated 

from  batch-to-batch.  Selection  of  batch  length  Q  represents 
a  compromise  between  fast  updates  to  counteract  interference 
non-stationarity  and  slow  updates  to  increase  secondary  data 
sample  support.  Recall  that  the  condition  KQ  >  2N  is 
required  to  provide  satisfactory  performance  for  the  SMI 
technique  [24],  The  adaptive  weights  wm  are  derived  from 
the  associated  integrated  sample  covariance  matrices  Rm, 

1  Qm— 1 

Rm  =  o  ^  (6) 

t=Q(m—l) 

and  applied  to  beamform  the  test  range  cell  in  the  current 
batch:  y(t)  =  w^x(f)  for  t  =  Q{m  —  1 ),..., Qm  —  1. 
The  similarity  between  the  matrices  R(f)  being  integrated 
depends  on  the  batch  length  Q  and  the  degree  of  interference 
spatial  non-stationarity.  In  general,  it  is  found  that  interference 
wavefronts  reflected  by  the  ionosphere  fluctuate  in  a  correlated 
manner  with  respect  to  time  [9],  [12],  In  simple  terms,  this 
means  that  the  matrices  R(t)  are  likely  to  be  similar  over  a 
number  of  adjacent  PRI,  but  become  progressively  different 
as  the  structure  of  the  interference  wavefronts  evolve  in  time 
(possibly  varying  appreciably  in  form  between  the  endpoints 
of  the  CPI).  A  parametric  interference  model  that  quantifies 
this  space-time  behaviour  was  experimentally  verified  in  [9], 
[13],  [18],  [19]  using  different  data  sets. 

B.  Surface  Clutter 

The  dominant  sea  clutter  contribution  is  produced  by  first- 
order  scatter  from  specific  spectral  components  of  the  ocean 
wavefield.  These  surface-height  components  are  called  Bragg 
waves.  Their  wavelength  A  is  exactly  half  the  radio  wavelength 
A  and  they  move  directly  towards  and  away  from  the  radar. 
In  water  of  depth  h  >  A/2,  the  Bragg  wave  trains  move 
with  radial  velocity  vg  =  ±^/gA/2Tr,  where  A  =  A/ 2  and 
g  is  acceleration  due  to  gravity.  This  imposes  a  Doppler  shift 
of  fb  =  ±s/ g/Xn  on  the  first-order  clutter  backscattered  at 
grazing  incidence  [20].  When  the  Bragg  wave  trains  propagate 
through  a  surface  current  with  mean  radial  velocity  vs,  an 
additional  Doppler  shift  fs  =  2vs/\  results  for  both  the 
advance  and  recede  Bragg  wave  spectral  returns. 

Let  nit )  be  the  first-order  ocean  clutter  received  at  the 
reference  antenna  element  in  the  array  due  to  a  scattering 
"patch”  defined  by  the  size  of  the  (unnamed)  test  range  cell 
and  a  narrow-beam  at  azimuth  9i.  In  accordance  with  [21], 
this  contribution  to  the  total  clutter  return  can  be  modelled  in 
the  time-domain  as 

r^t)  =  ^4ie-t427r(/s+/b)t+<)>i(t)}  +  A.ej{2n(fa-fb)t+<l>i(t)}  ^ 

where  At  and  A ,  are  the  advance  and  recede  first-order  clutter 
amplitudes  respectively.  These  amplitudes  are  proportional  to 
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the  ocean  directional  wave-height  spectrum  at  the  Bragg  wave- 
vectors  (generally  two  orders  of  magnitude  higher  than  the 
surrounding  clutter  continuum  due  to  higher  order  scatter). 
As  the  primary  concern  is  not  to  destroy  the  natural  temporal 
coherence  of  the  first-order  clutter,  we  focus  on  modelling 
these  dominant  returns  and  omit  the  much  weaker  higher 
order  scatter.  In  practice,  the  surface-current  through  the  radar 
footprint  may  have  a  non-uniform  velocity  component  with 
different  dynamic  behaviour  from  one  scattering  patch  to 
another.  Surface-current  turbulence  over  the  data  collection 
period  imposes  a  phase  modulation  <j>i(t)  on  the  Bragg  wave 
spectral  components,  a  mechanism  that  broadens  the  approach 
and  recede  Bragg  lines  into  identical  shapes  [22], 

The  overall  first-order  ocean  clutter  received  at  different 
antenna  elements  is  the  vector  addition  of  the  narrow-beam 
returns  from  I  azimuth  cells  subtending  the  intersection  of  the 
transmit  and  receive  beampatterns:  c(t)  =  ri{t)'v{6i).  If 

a  deterministic  ground  clutter  component  cg  with  zero  Hertz 
Doppler  shift  is  also  included,  ignoring  the  movement  of 
ground  scatters  and  non-ideal  systems  effects,  it  follows  that 
the  clutter  snapshots  received  by  the  array  can  be  modelled  as 

c(t)  =  cg  +  ca(ty2”(f‘+fb)t  +  c  r(t)ej2w(-fs~fb)t  (8) 

where  c a(t)  =  J2i=o  Aje-^^Wv^)  is  the  modulation  on 
the  advance  wave  and  c r(t)  =  J2i=o  ^4iej27r?ii^v(#i)  is  that 
for  the  recede  wave.  The  instantaneous  Bragg  frequencies 
change  relatively  slowly  and  may  be  regarded  almost  constant 
over  a  sufficiently  short  sub-CPI.  The  assumption  of  “frozen” 
Bragg  frequencies  over  Q  consecutive  PRI  allows  the  clutter 
to  be  represented  as  superposition  of  complex  sinusoids  with 
fixed  frequencies  over  each  sub-CPI.  Such  a  model  has  previ¬ 
ously  been  adopted  to  parametrically  estimate  surface  currents 
[21]  and  to  cancel  clutter  over  CPI  lengths  shorter  than  3 
seconds  [23]. 

In  mathematical  terms,  we  let  ca(t)  =  ca(TTO)ej27r‘5^m<  and 
cr(t)  =  cr(rm)ej27r,5^mt  for  t  =  Q(m—  1), . . . ,  Qm  —  1  where 
Tm  =  Q(m  —  1)  and  6fm  is  the  slowly  changing  shift  in  Bragg 
wave  instantaneous  frequency  that  causes  incoherent  integra¬ 
tion  of  the  first-order  ocean  clutter  over  the  CPI.  By  defining 
the  matrix  Am  =  [cs,  ca(rm),  cr(rm)]  and  the  vector  p(t)  = 
llte^2n^m+^t,e^2n^m~^t]T  for  fm  =  fs  +  Sfm,  the  first 
order  multi-channel  clutter  realisation  is  locally  represented  by 
a  dynamic  spatial  subspace  model  with  an  instantaneous  rank 
of  L  =  3  over  each  sub-CPI. 

c(f)  =  Amp(f)  ,  t  =  Q(m  -  1), ... ,  Qm  -  1  (9) 

In  general,  the  minimum  value  of  L  is  set  by  the  number  of 
dominant  components  resolved  in  the  clutter  Doppler  spectra. 
For  example,  L  =  3  corresponds  to  two  resolved  Bragg  lines 
and  ground  clutter.  Values  of  L  >  3  may  be  adopted  to  account 
for  subspace  leakage  caused  by  the  continuous  nature  of  the 
actual  phase  modulation  process  or  to  account  for  some  of  the 
second-order  clutter.  The  condition  L  <  Min(Q,  N )  must  be 
satisfied  since  Q  snapshots  of  dimension  N  can  be  represented 
by  a  subspace  of  dimension  less  than  or  equal  to  Min(Q,  N ). 


III.  Proposed  Algorithm 

The  first  adaptive  beamformer  wj  in  the  sequence  is  chosen 
to  provide  maximum  SINR  for  the  integrated  sample  covari¬ 
ance  matrix  Ri  in  accordance  with  the  minimum  variance 
distortionless  response  (MVDR)  criterion.  This  weight  vector 
is  given  by  wi  =  [v**  (0)Rj"1v(0)]_1Rj"1v(0)  and  used  to 
beamform  the  array  snapshots  x(/:j  received  in  the  first  batch 
at  the  test  range  cell;  y(t)  =  w^x(f)  for  t  =  0, . . . ,  Q  —  1. 

To  avoid  the  coherence  “discontinuity”  across  the  batches 
without  prohibiting  the  adaptive  weight  vector  to  change  from 
batch-to-batch,  the  second  weight  vector  W2  is  formed  in 
a  manner  that  is  dependent  on  the  first  wj  such  that  the 
instantaneous  change  in  antenna  pattern  is  constrained  to  be 
orthogonal  to  the  clutter  subspace  in  the  second  batch.  Stated 
mathematically,  the  key  idea  of  the  proposed  algorithm  is  to 
closely  approximate  the  condition  (w2  —  W|  )H  A.2  =  0,  as 
this  permits  W2  to  change  from  wj  for  interference  rejection 
purposes  while  at  the  same  time  regulating  this  change  such 
that  it  is  not  experienced  by  the  clutter  c(t')  =  A2p(f)  received 
in  the  second  batch. 

The  clutter  subspace  in  each  batch  is  not  observable  because 
x(t)  contains  interference  and  noise.  However,  the  adaptive 
filter  wi  is  still  effective  for  interference  removal  in  the 
first  L  pulses  of  the  second  batch.  If  necessary,  this  can 
be  ensured  by  extending  the  covariance  matrix  estimate  of 
Ri  =  (Q  +  L)_1  1  R-M  t0  include  the  first  L  pulses 

of  second  batch.  Applying  this  filter  to  the  first  L  snapshots 
of  the  second  batch  yields  wfrx(f)  ss  g(t)  +  w frc(f)  for 
t  =  Q, . . .  ,Q  +  L  —  1,  where  the  residual  interference-plus- 
noise  contribution  at  the  beamformer  output  is  assumed  to 
be  negligible  compared  with  the  unrejected  clutter.  Using  the 
same  argument,  wf^x(i)  ss  <7(i)  +  w|*c(f)  for  t  =  Q, . . . ,  Q  + 
L  —  1  since  both  weight  vectors  have  unity  gain  response  to  the 
useful  signal.  If  the  second  adaptive  filter  is  made  to  satisfy  L 
linear  data-driven  constraints  defined  by  w^x(f)  =  w^x(f) 
for  t  =  Q, . . . ,  Q  +  L  —  1,  then  to  a  good  approximation,  the 
clutter  contribution  at  the  output  of  the  second  filter  is  as  if 
it  were  beamformed  by  the  first  filter  (i.e.  w|rc(f)  «  w^c(f) 
is  approximated  for  t  =  Q, . . .  ,Q  +  L  —  1)  even  though  the 
two  filters  may  be  quite  different  (wi  /  W2),  as  required  for 
effective  interference  rejection. 

As  the  L  clutter  snapshots  c(t)  =  A2p(f)  at  t  =  Q, . . . ,  Q+ 
L  —  1  are  linearly  independent,  they  collectively  span  the  col¬ 
umn  space  of  the  matrix  A2.  Consequently,  the  L  data-driven 
constraints  imposed  on  W2  encourage  the  relative  change  in 
antenna  pattern  to  be  orthogonal  to  the  clutter  subspace  (W2  — 
wi)hA2  ~  0.  This  condition  carries  with  it  the  implication 
that  the  clutter  contribution  at  the  output  of  the  second  filter  is 
as  if  it  were  beamformed  by  the  first  filter  over  the  whole  of  the 
second  batch  t  =  Q, . . . ,  20  1  despite  the  use  of  only  L  data- 
driven  constraints.  By  defining  C2  =  [v(0),x(Q), . . .  ,x((3  + 
L  -  1)]H  and  f2  =  [l,xff(Q)wi, . .  .,xH(Q  +  L  -  l)wi]T 
for  the  second  batch,  W2  is  synthesised  according  to  (10). 

minwfrR2W2  subject  to  C2W2  =  f2  (10) 

W2 
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This  is  the  linearly  constrained  minimum  variance  (LCMV) 
optimisation  problem  which  has  the  following  solution  [25]. 

w2  =  R^1C2[CfR^1C2]“1f2  (11) 

The  remaining  batches  comprising  the  CPI  are  processed  in 
similar  fashion  by  iterating  with  respect  to  m  =  2 , ,M 
such  that  the  ?nth  weight  vector  wm  is  computed  as, 

wm  =  R”1Cm[C"R”1Cm]-1fm  (12) 

for  Cm  =  [v(0),x(Q(m  -  1)), . . .  ,x(Q(m  -  1)  +  L  -  1)}H 
and  fm  =  [1  ,xH(Q(m-  l))wm_i, . . . ,  xH {Q(m  -  1)  +  L  - 
l)wm_i]T.  That  is,  the  iterations  are  continued  until  the  whole 
CPI  is  processed.  Note  that  the  proposed  algorithm  is  iterative 
in  the  sense  that  the  MVDR  weight  vector  wj  initialises  the 
sequence  at  m  =  1  and  the  weights  for  batches  m  =  2, . . . ,  M 
are  computed  using  the  known  weights  from  the  previous  step. 

The  adopted  first-order  clutter  model  allows  the  proposed 
algorithm  to  update  the  weight  vector  in  non-overlapping 
batches  whose  length  is  primarily  determined  by  the  degree 
of  interference  spatial  non-stationarity.  This  is  different  to 
the  SC  approach,  where  the  adaptive  filter  is  re-adjusted 
in  “sliding  window’’  fashion  every  PRI  (irrespective  of  the 
interference  characteristics)  in  an  attempt  to  protect  the  AR 
clutter  Doppler  spectrum  properties.  The  new  philosophy  not 
only  provides  robustness  against  the  accumulation  of  weight 
vector  estimation  errors  over  the  CPI,  but  also  dramatically 
reduces  the  number  of  matrix  inversions  required  and  hence 
the  computation  load.  For  example,  if  P  =  256,  L  =  3 
and  M  =  16,  the  proposed  method  reduces  the  number 
of  N  x  N  sample  covariance  matrix  inversions  by  a  factor 
F  =  pfIL  «  16  which  yields  appreciable  computational 
benefits  for  enabling  real-time  implementation. 

IV.  Experimental  Results 

Experimental  data  for  this  study  was  collected  using  the 
Iluka  HFSW  radar,  located  near  Darwin  in  far  north  tropical 
Australia.  The  radar  is  bistatic  with  a  high  power  (1-10  kW) 
transmit  site  at  Stingray  Head  (65  km  south-west  of  Darwin) 
and  a  lower  power  (100  W)  site  at  Lee  Point  (10  km  north¬ 
east  of  Darwin).  The  receiving  system  at  Gunn  Point  (30  km 
North-East  of  Darwin)  is  based  on  a  500  m  long  ULA  of  32 
vertical  monopoles,  each  antenna  element  is  connected  to  a 
HF  receiver  and  two  dummy  elements  are  added  at  either  end 
of  the  array  to  reduce  the  effects  of  mutual  coupling. 

Each  CPI  is  32  seconds  long  and  consists  of  P  =  256  linear 
frequency  modulated  continuous  waveform  (FMCW)  pulses 
or  sweeps  emitted  with  centre  frequency  fc  =  7.719  MHz, 
bandwidth  B  =  50  kHz  and  pulse  repetition  frequency  (PRF) 
fp  =  8  Hz.  The  analysis  is  based  on  a  ULA  of  N  20  well- 
calibrated  receivers,  the  full  array  could  not  be  used  because 
of  some  poor  receivers  that  produced  unreliable  output.  A 
total  of  60  range  cells  were  retained  after  receiver  mixing 
and  filtering,  the  first  12  clutter- free  range  cells  providing  the 
secondary  snapshots  for  training  the  adaptive  beamformers. 
The  interference  background  was  unknown  and  possibly  arose 
due  to  a  multiplicity  of  man-made  and  natural  sources. 


A.  Conventional  and  SMI-MVDR  Beamforming 

Using  Q  =  16  PRI  and  M  =  16  batches  (i.e.  update  interval 
of  2  seconds),  the  SMI-MVDR  beamformer  was  formed  using 
I\  =  10  range  cells  per  PRI  (KQ  =  160  =  8 N)  and  applied 
to  clutter-free  range  cell  12  to  investigate  interference  rejection 
across  ranges.  Figure  1  compares  the  interference-plus-noise 
Doppler  spectra  processed  by  the  conventional  beamfomer 
(Hamming  taper)  and  the  time-varying  SMI-MVDR  adaptive 
beamformer  with  no  data-driven  constraints.  Although  both 
beamformers  are  normalised  to  unit  gain  in  the  look  direction, 
the  time-varying  adaptive  beamformer  is  seen  to  improve 
interference  rejection  by  15-20  dB  across  Doppler  space. 

Figure  2  shows  the  application  of  the  same  weights  to 
beamform  range  cell  16  which  additionally  contains  clutter 
and  a  real  target  (aircraft  echo)  in  the  vicinity  of  Doppler 
bin  38.  Although  the  time-varying  SMI-MVDR  beamformer 
effectively  rejects  interference,  it  severely  degrades  target 
detection  because  the  changing  antenna  patterns  destroy  the 
temporal  coherence  of  the  clutter  (particularly  the  first-order 
returns)  which  smear  across  Doppler  space  and  totally  obscure 
the  target  echo. 

B.  Proposed  Algorithm  with  Data-Driven  Constraints 

The  same  data  was  processed  by  the  proposed  adaptive 
beamformer  using  L  =  1,2  data-driven  constraints.  Figure  3 
shows  that  L  —  1  performs  very  poorly,  a  choice  of  L  =  2 
makes  the  target  slightly  more  distinguishable  but  severely 
disturbs  the  clutter  Doppler  spectrum  compared  with  the 
conventional  beamformer  output.  It  is  evident  from  Figure  4 
that  the  use  of  L  =  3  data-driven  constraints  improves  the 
SINR  of  a  real  target  by  approximately  20  dB  relative  to  the 
conventional  beamformer  with  no  noticeable  degradation  in 
sub-clutter  visibility.  In  other  words,  this  time-sequence  of 
adaptive  weight  vectors  is  felt  by  the  interference  but  remains 
effectively  “invisible”  to  the  clutter. 

Three  ideal  synthetic  targets  were  injected  in  Doppler 
bins  53,  112  and  218  at  test  range  cell  33  to  compare  the 
proposed  adaptive  beamformer  with  the  steady-state  or  “fixed” 
SMI-MVDR  beamformer  trained  on  the  same  I\  =  10  range 
cells  but  over  the  whole  CPI  (i.e.  Q  =  256  PRI).  Figure  5 
shows  that  while  both  approaches  can  detect  strong  signals, 
either  in  the  noise  (bin  53)  or  near  the  clutter  (bin  112),  only 
the  proposed  method  is  capable  of  detecting  the  weaker  high 
velocity  target  (bin  218). 

C.  Comparison  with  Stochastic  Constraints  Method 

Figure  6  compares  the  performance  of  the  proposed  adaptive 
beamformer  with  L  =  3  data-driven  constraints  against  the 
SC  method  [  14]  based  on  a  third  order  AR  clutter  model.  The 
former  is  robust  and  produces  sharper  Bragg  lines  that  enable 
a  low  velocity  target  near  the  right  Bragg  line  (bin  112)  to 
be  detected.  The  large  number  of  updates  in  the  SC  method 
results  in  a  significant  accumulation  of  estimation  errors  that 
not  only  widens  the  Bragg  lines,  but  also  smears  clutter  energy 
into  other  regions  of  Doppler  space,  making  the  weak  high 
velocity  target  (bin  218)  more  difficult  to  detect. 
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Apart  from  the  relative  improvement  in  performance,  the 
proposed  adaptive  beamfomer  is  about  three  times  faster  to 
compute  than  the  SC  method  (timing  tests  run  on  an  Alpha 
(DS20)  500  MHz  computer).  Hence,  the  proposed  algorithm 
represents  an  attractive  candidate  for  real-time  implementation 
in  practical  HFSW  radar  systems. 

V.  Summary  and  Conclusion 
This  paper  has  presented  a  robust  and  computationally 
efficient  adaptive  beamforming  algorithm  for  the  cancellation 
of  spatially  non- stationary  interference  while  preserving  the 
Doppler  spectrum  characteristics  of  the  dominant  (first-order) 
clutter  returns.  The  new  method  updates  the  weight  vector  at  a 
rate  determined  the  degree  of  interference  non-stationarity,  this 
is  in  contrast  with  the  SC  approach  which  re-adjusts  the  weight 
vector  every  PRI  regardless  of  the  interference  characteristics 
in  an  attempt  to  protect  the  AR  Doppler  spectrum  properties 
of  the  clutter  (an  update  rate  that  is  often  faster  than  necessary 
for  effective  HF  interference  rejection). 

Experimental  results  confirm  that  the  proposed  method  can 
improve  the  SINR  of  real  targets  by  20  dB  relative  to  the 
conventional  beamformer  in  a  HFSW  radar.  Apart  from  the 
computational  advantage  with  respect  to  the  SC  technique, 
the  proposed  method  was  found  to  be  more  robust  against  the 
accumulation  of  estimation  errors  which  degrade  sub-clutter 
visibility.  The  benefits  of  the  new  approach  with  respect  to 
previously  existing  methods  have  been  demonstrated  on  a 
HFSW  radar,  although  it  is  envisaged  that  the  same  technique 
can  be  applied  to  skywave  OTH  and  other  types  of  radar. 
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Fig.  1.  Interference-plus-noise  Doppler  spectra  for  conventional  (Hamming  Fig.  4.  Operational  radar  Doppler  spectra  for  the  conventional  (Hamming 
taper)  beamformer  and  time-varying  MVDR  adaptive  beamformer  with  no  taper)  beamformer  and  proposed  adaptive  beamformer  with  L  =  3  data-driven 
data-driven  constraints.  constraints. 
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Fig.  2.  Operational  radar  Doppler  spectra  for  the  conventional  (Hamming  Rg  5  Operational  radar  Doppler  spectra  with  synthetically  injected  targets 

taper)  beamformer  and  time-varying  MVDR  adaptive  beamformer  with  no  for  the  fixed  MVDR  and  propoSed  adaptive  beamformers. 

data-driven  constraints. 
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Fig.  3.  Operational  radar  Doppler  spectra  for  the  proposed  adaptive 
beamformer  with  L  =  1  and  L  —  2  data-driven  constraints. 
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Fig.  6.  Operational  radar  Doppler  spectra  with  synthetically  injected  targets 
for  the  stochastic  constraints  and  proposed  adaptive  beamformers. 
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